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Energies of {3s'^3p'^3(fmi'), {3s'^3p^3d^°mi'), and (3s3p^3d^°4M/') states for Cu-like ions with 
Z = 30 -100 are evaluated to second order in relativistic many-body perturbation theory (RMBPT) 
starting from a Ni-like Dirac-Fock potential. Second-order Coulomb and Breit-Coulomb interactions 
are included. Correction for the frequency-dependence of the Breit interaction is taken into account 
in lowest order. The Lamb shift correction to energies is also included in lowest order. Intrinsic 
particle-particle-hole contributions to energies are found to be 20-30% of the sum of one- and two- 
body contributions. Transition rates and line strengths are calculated for the 31 — 41' electric-dipole 
(El) transitions in Cu-like ions with nuclear charge Z = 30 - 100. RMBPT including the Breit 
interaction is used to evaluate retarded El matrix elements in length and velocity forms. First-order 
RMBPT is used to obtain intermediate coupling coefficients and second-order RMBPT is used to 
calculate transition matrix elements. A detailed discussion of the various contributions to the dipole 
matrix elements and energy levels is given for copperlike tungsten {Z = 74). The transition energies 
used in the calculation of oscillator strengths and transition rates are from second-order RMBPT. 
Trends of the transition rates as functions of Z are illustrated graphically for selected transitions. 
Comparisons are made with available experimental data. These atomic data are important in 
modeling of M-shell radiation spectra of heavy ions generated in electron beam ion trap experiments 
and in M-shell diagnostics of plasmas. 

PACS numbers; 31.15.Ar, 31.15.Md, 31.25.Jf, 32.30.Rj 



I. INTRODUCTION 

This work further develops the apphcation of the rela- 
tivistic many-body perturbation theory (RMBPT) to the 
studies of atomic characteristics of particle-particle-hole 
excitations of closed-shell ions. Recently, RMBPT calcu- 
lations of energies ,1] and transition rates in Na-like 
ions have been performed. The present paper focuses on 
the RMBPT calculations of energies and transition rates 
in Cu-like ions. These transitions form satellite lines to 
the brightest Ni-like ions and are of a great importance 
for M-shell diagnostics of heavy ions. 

The second-order RMBPT calculations for Cu-likc ions 
start from a 1522522^63^23^63^10 [js^jj Dirac-Fock poten- 
tial. All possible 31 holes and AlAl' particles leading to 
563 odd-parity and 571 even-parity 3Z~^4r4r'(J) states 
are considered. The energies of the 1134 states and tran- 
sition probabilities of the 2294 electric-dipole lines are 
calculated for Cu-like ions with Z ~ 30-100. Transi- 
tions from the (Ss^Sp^Scfmi'), {Ss^Sp^Sd^^mi'), and 
(353/3^104/4/') states to the ground (3523/3^1043) qj. 
singly excited {Ss'^Sp^Sd^^il) states (with / — p,d, f) 
form satellite lines to the bright electric-dipole (El) lines 
in Ni-like ions. These core-excited states (or often called 
doubly-excited states) in copperlike ions have been stud- 
ied extensively both experimentally and theoretically in 
the past 15-20 years. 

Experimentally, these configurations have been stud- 
ied by photon and electron emission spectroscopy. To 



our knowledge, the first measurements of 3d — Ap tran- 
sitions in Cu-like W and Tm were done by Klapisch 
et al. by classification of x-ray spectra from laser 
produced plasmas in the range 6 - 9 A. It was shown 
that most of the Cu-like 3d — 4p line radiation came 
from (3s23/3di°4s - 3s^3p^3d^As'ip), {3s^3p^3d^^Ap - 
3523/3^^4/), (3s23/3di04d - 3s23/3d94p4d), and 
(3523/3(ii°4/ - 3523/3^94^4/) transitions. Wave- 
lengths and transition probabilities were calculated in Q 
by the relativistic parametric potential method . The 
same method was extended by Mandelbaum et al. 
to study x-ray spectra from laser produced plasmas of 
atoms from Tm (Z = 69) up to Pt {Z = 78). An ex- 
tended analysis of x-ray spectra of laser produced gold 
plasmas has been performed by Busquet et al. Q in- 
cluding 3(i — 4/, 3d — Ap, and 3p — As transitions in Ni-, 
Co-, Cu-, Zn-, and Ga-like Au. The Ni-, Co-, Cu-, Zn-, 
Ga-, Ge-, and As-like isoelectronic sequences have been 
considered in Ref. to investigate the x-ray spectra of 
laser-irradiated elements from tantalum {Z = 73) to lead 
{Z = 82). In addition to the above mentioned isoelec- 
tronic sequences, the Fe- and Mn-like states have been 
included by Zigler et al. 0| to analyze x-ray spectra em- 
mited by laser-produced plasmas of lanthanum {Z = 57) 
and praseodymium {Z = 59). The wavelengths and tran- 
sition probabilities have been calculated in 0,0, Q by 
the relativistic parametric potential method |4|. Inves- 
tigation of the x-ray spectrum emitted by laser produced 
bariurnplasmas has been recently described by Doron 
et al. [a US- The RELAC relativistic atomic code [J 
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has been used to identify 3d — nl (n = 4 to 8), 3p — 4s, 
and 3p — 4d transitions of Ni-like Ba^*+ and correspond- 
ing satellite transitions in neighboring ionization states: 
Fe-, Co-, Cu-, Zn-, Ga-, and Ge-like ions. The X-ray 
spectrum emitted by a laser produced cerium plasma in 
the range 7.5 - 12 A has been recently investigated in 
detail by Doron et al. [HIS- The RELAC computer 
code Q has been used to study x-ray spectra from highly 
charged t ung sten ions in tokamak plasmas in the range 
7 - 10 A [13. 

In the present paper, RMBPT is implemented to de- 
termine energies of {3s^3p^3(fmi'), {3s^3p^3d^°mi'), 
and {3s3p^3d^'~' AlAV) states for Cu-like ions with nuclear 
charges in the range of Z = 30 -100. The calculations are 
carried out to second order in perturbation theory and in- 
clude second-order Coulomb and Breit interactions. Cor- 
rections for the frequency-dependent Breit interaction are 
taken into account in the lowest order. Screened self- 
energy and vacuum polarization data given by Kim et al. 
[T^ are used to determine the QED correction. 

RMBPT is used to evaluate matrix elements, line 
strengths, and transition rates for 2294 allowed and for- 
bidden electric-dipole transitions between the 1332 even- 
parity core-excited states [Ni](3d~^ -|- 3s~^)(4s^ -I- 4s4d-|- 
ip'^ + M'^ + ApAf+Ap) + [Ni](3p-i)(4s4p-|-4s4/-|-4p4d-h 
4(i4/), and the singly excited [Ni](4p-|-4/) states and the 
962 odd-parity core-excited states [Ni](3p~^)(4s2-|-4s4(i-t- 
Ap^ + Ad^ + ApAf + ip) + [Ni](3d-i+3s-i)(4s4p+4.s4/+ 
ApAd + AdAf) and the ground state [Ni](4s) together with 
the singly excited [Ni] (Ad) states in Cu-likc ions with nu- 
clear charges ranging from Z = 30 to 100. Retarded El 
matrix elements are evaluated in both length and veloc- 
ity forms. These calculations start from a [Ni] Dirac- 
Fock potential. First-order perturbation theory is used 
to obtain intermediate coupling coefficients and second- 
order RMBPT is used to determine transition matrix el- 
ements. The transition energies employed in the calcu- 
lations of line strengths and transition rates are derived 
from second-order RMBPT. 



II. METHOD 

Details of the RMBPT method were presented in 
Ref. for calculation of energies of particle-particle-hole 
states and in Ref. '5| for calculation of radiative transi- 
tion rates from particle-particle-hole state to one-particle 
state. Differences between calculations for Na-like and 
Cu-like ions are due to the increased size of the model 
space {Al'Al"3l-^ instead of 3l'3l"2l-^) and differences 
in the Dirac-Fock potential (1522522^63323^63^10 instead 
of ls^2s^2p^), leading to 1134 states instead of 106 and 
more laborious numerical calculations. 

As a first step, we determine and store the single- 
particle contributions to the energies for five n=3 hole 
states (3s, 3pi/2, 3p3/2, 3^3/2, and 3^5/2) and the 
seven n=A valence states (45, 4pi/2j 4p3/2, 4^3/2, 
4(^5/2 J 4/5/2, and 4/7/2) in lowest, first, and second or- 



ders. Next, we evaluate and store the 664 two-particle 
{AlAl'J\H''^\Al"Arj) matrix elements and the 1127 hole- 
particle {31AI' J\H''^\3l"Al"' J) matrix elements of the ef- 
fective Hamiltonian in first and second orders. It should 
be noted that these one-particle, two-particle, and hole- 
particle matrix elements were used previously to evaluate 
energies of the AlAl' levels in Zn-like ions 0| and energies 
of the 3l~^Al' levels in Ni-like ions |^^. Finally, second- 
order particle-particle-hole matrix elements are evaluated 
p. Combining these data using the method described 
below, we calculate one-, two-, and three-body contribu- 
tions to the energies of Cu-like ions. 

The calculations are carried out using sets of basis 
Dirac-Hartree-Fock (DHF) orbitals. The orbitals used 
in the present calculation are obtained as linear combi- 
nations of B-splines. These B-spline basis orbitals are de- 
termined using the method described in Ref. 0|. Forty 
B-splines of order of eight for each single-particle angu- 
lar momentum state are used and all orbitals with orbital 
angular momentum I < 7 are included in the basis set. 



A. Model space 

The model space for core-excited AlAl'3l~^ states of 
copperlike ions includes 563 odd-parity states consisting 
of 78 J=l/2 states, 131 J=3/2 states, 143 J=5/2 states, 
125 J=7/2 states, and 86 J=9/2 states. Additionally, 
there are 571 even-parity states consisting of 78 J =1/2 
states, 131 J=3/2 states, 148 J=5/2 states, 125 J=7/2 
states, and 89 J=9/2 states. The distribution of some 
of the 1134 states in the model space is summarized in 
Table I of the accompanying EPAPS document p^. 



B. Energy-matrix elements 

The evaluation of the second-order energies for the 
AlAl'{Ji)3l" (J) states in Cu-like ions follows the pat- 
tern of the corresponding calculation for Zn-like and Ni- 
like ions given in Refs. J^S^, ,16]. In particular, we use 
the second-order one- and two-particle matrix elements 
for Zn-like ions calculated in [13 and hole-particle ma- 
trix elements for Ni-like ions calculated in [l^, but re- 
coupled as described below, to obtain the one- and two- 
particle contributions for Cu-like ions. We will discuss 
how these matrix elements are combined to obtain the 
one- and two-particle contributions to energies of Cu-like 
ions. We refer the reader to Ref. fislll^ for a discussion 
of the how the basic one- and two-particle matrix ele- 
ments were evaluated. An intrinsic particle-particle-hole 
diagram also contributes to the second-order energy for 
Cu-like ions. It should be noted that the angular part 
of the expression for the particle-particle-hole diagram 
differs from the corresponding expression for the three- 
particle matrix elements given in Ref. ^ detailed 
discussion of this difference is given in Ref. [ll]. 
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TABLE L Second-order contributions to tiie energy matrices (a.u.) for odd-parity states with J=l/2 in the case of Cu-like 
tungsten, Z=74. One-body, two-body, ant three-body second-order Coulomb contributions are given in columns labelled E\ , 
and eI^\ respectively. 



4/ljl4Z2j2[Jl2]3;3 


4/lil4Z2j2[Jl2]3;3 


El ' 


n(2) 
E2 


n(2) 
-^3 


n(2) 


4Sl/24pi/2(l)3d3/2 


4Sl/24pi/2(l)3d3/2 


-0.247880 


0.104563 


0.025088 


-0.118235 


4Sl/24p3/2 (2)3^5/2 


4Sl/24p3/2(2)3d5/2 


-0.226244 


0.113517 


0.020801 


-0.091926 


4Sl/24p3/2(l)3d3/2 


4Sl/24p3/2(l)3d3/2 


-0.237415 


0.076978 


0.029703 


-0.130734 


4Sl/24p3/2(2)3d3/2 


4Sl/24p3/2(2)3d3/2 


-0.237415 


0.124305 


0.021436 


-0.091674 


4pi/24d3/2(2)3d5/2 


4pi/24d3/2(2)3d5/2 


-0.240289 


0.122424 


0.028030 


-0.089836 


4pi/24d5/2(2)3d5/2 


4pi/24d5/2(2)3d5/2 


-0.236161 


0.031594 


0.009588 


-0.194980 


4pi/24d5/2(3)3d5/2 


4pi/24d5/2(3)3d5/2 


-0.236161 


0.001311 


-0.002282 


-0.237132 


4pi/24d3/2 (1)3^3/2 


4pi/24d3/2(l)3d3/2 


-0.251460 


0.040039 


0.036490 


-0.174931 


4pi/24d3/2(2)3d3/2 


4pi/24d3/2(2)3d3/2 


-0.251460 


0.060163 


0.026395 


-0.164902 


4Sl/24/5/2(2)3d5/2 


4Sl/24/5/2(2)3d5/2 


-0.229530 


0.067699 


0.039638 


-0.122194 



TABLE IL Energies of selected odd-parity levels with J=l/2 of Cu-like tungsten, Z=74 in a.u. = Eq + Ei + Bi. 



jj coupling 




Bi 


E2 


-Elamb 


-Etot 


4sl/24pi/2(l)3d3/2 


-25.954298 


-0.013561 


-0.112483 


0.069035 


-26.011307 


4Sl/24p3/2(2)3d5/2 


-24.879634 


-0.005093 


-0.090398 


0.065104 


-24.910021 


4Sl/24p3/2(l)3d3/2 


-22.639268 


-0.054704 


-0.090313 


0.073682 


-22.710603 


4Si/24p3/2(2)3d3/2 


-21.914637 


-0.057899 


-0.125159 


0.073518 


-22.024177 


4pi/24d3/2(2)3d5/2 


-16.127594 


0.032360 


-0.088446 


-0.001127 


-16.184807 


4pi/24d5/2(2)3d5/2 


-14.892864 


-0.001171 


-0.108506 


0.001320 


-15.001221 


4pi/24d5/2(3)3d5/2 


-14.149853 


0.000699 


-0.136464 


0.003373 


-14.282246 


4pi/24d3/2(l)3d3/2 


-13.066261 


-0.033638 


-0.123168 


0.007346 


-13.215720 


4pi/24d3/2(2)3d3/2 


-12.799288 


-0.034773 


-0.108188 


0.050678 


-12.891570 


4Sl/24/5/2(2)3d5/2 


-12.654131 


-0.048963 


-0.119286 


0.050671 


-12.771709 



Table (see also Tabic II of Ts^ ) provides an illustra- 
tion of various contributions to the second-order energies 
for the special case of Cu-like tungsten, Z — 74. In this 
table, we show the one-body, two-body and three-body 
second-order Coulomb contributions to the energy matrix 

labelled as ^ » = 1, 2, 3, The one-body second-order 
(2) 

energy, E{ is obtained as the sum of the two valence and 

one hole Ei^^ energies. The values of ' are non-zero 
only for diagonal matrix elements. Even for odd-parity 
states with J =1/2 there are 78 diagonal and 6006 non- 
diagonal matrix elements for AlAl'{Ji)'il" (J) core-excited 
states. We list data only for the first ten diagonal matrix 
elements of odd-parity states with J=l/2 in Tabled (a 
more complete set of data is given in Table II of J^] ) . It 
can be seen from the table that two-body and three-body 
second-order contributions are positive, when the one- 
body contributions are negative. The three-body contri- 
butions give about 20% in the total second-order contri- 
butions. The values of the £"3 and E2 non-diagonal 
contributions are smaller than values of diagonal contri- 
butions by factor of 3-5. 

After evaluating the energy matrices, eigenvalues and 
eigenvectors are calculated for states with given values 
of J and parity. There are two possible methods to 
carry out the diagonalization: (a) diagonalize the sum 
of zeroth- and first-order matrices, then calculate the 



second-order contributions using the resulting eigenvec- 
tors; or (b) diagonalize the sum of the zeroth-, first- and 
second-order matrices together. Following Ref. 0, we 
choose the second method here. 

Tablelnllists the energies of ten excited states of W^^+ 
from Tabled including the total energies i?tot- The latter 
is the sum of the following contributions: = E'^^^ 

+E'''^^ +B^^\ the second-order Coulomb and Breit energy 
and the QED correction Slamb- The QED correc- 
tion is approximated as the sum of the one-electron self 
energy and the first-order vacuum-polarization energy. 
The screened self-energy and vacuum polarization data 
given by Kim et al. M] are used to determine the QED 
correction i?LAMB (see, for detail Ref. 0). As can be 
seen, the second order correction contributes to the total 
energy from 0.4 % for lowest levels up to 8 % for high- 
excited levels. The levels in this table (see also Table 
III of 0) could be divided into groups corresponding to 
excited 4141' states and 31 hole states. 



C. Dipole matrix element 

We designate the first-order dipole matrix element by 
Z^-^\ the Coulomb-correction to the second-order matrix 
element Z^^\ and the second-order Breit correction B'-^K 
The evaluation of Z^-^\ Z^^\ and B'-^^ for Cu-Hke ions 
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TABLE III: Uncoupled reduced matrix elements in length L and velocity V forms for transitions between the selected odd-parity 
core-excited states with J — 1/2 and the ground 4s and singly excited 4d3/2 states in W*^^ ion. 



AljAl J (Ji)3/ J 


7(1) 


7(1) 




Zv 




c(2) 


Pl 


p(dcrv) 






[mi'nJiW'j" 


(1/2)- 4si/2] 


transitions 








4Sl/24pi/2(l)3d3/2 


0.047738 


0.044870 


0.002168 


-0.000123 


-0.000050 


0.001693 


0.047424 


0.000014 


4Sl/24p3/2(2)3d5/2 


0.052275 


0.049040 


0.003284 


0.000010 


-0.000034 


0.002425 


0.051930 


-0.000079 


4Sl/24/5/2(2)3d5/2 


0.041125 


0.038959 


-0.001778 


0.000268 


0.000087 


-0.000826 


0.041123 


0.000270 


4Sl/24/5/2 (3)3^5/2 


0.034757 


0.032927 


-0.001517 


0.000227 


0.000074 


-0.000842 


0.034756 


0.000225 


4Sl/24d3/2(l)3pi/2 


-0.067725 


-0.064364 


0.000594 


0.000439 


-0.000594 


0.000216 


-0.066887 


-0.000685 


^•'l/2^P3/2l-'-j'J«l/2 


n4'^7'i1 


-0.041541 


-0.017203 








n4984'i 

U.U4:^04:tJ 








[mi'nj^w'j" 


(1/2)- 4d3/2] 


transitions 








4pi/24d3/2(2)3d3/2 


0.046223 


0.043446 


0.002084 


-0.000119 


-0.000054 


0.001629 


0.045918 


0.000010 


4p3/24d3/2(3)3d5/2 


0.048898 


0.045873 


0.002691 


0.000010 


-0.000028 


0.002093 


0.048576 


-0.000065 


4Sl/24(i3/2(2)3p3/2 


0.062272 


0.058739 


0.001399 


-0.000036 


0.000201 


0.001391 


0.061840 


0.000223 


4d3/24/7/2 (2)34/2 


0.189383 


0.179458 


-0.008551 


-0.001590 


0.000298 


-0.003992 


0.187961 


0.000470 


4^3/24/7/2 (3)3d5/2 


-0.146695 


-0.139008 


0.006235 


0.001231 


-0.000236 


0.003225 


-0.145594 


-0.000376 


4^3/24/5/2(1)3^3/2 


0.146008 


0.138475 


-0.007222 


-0.000623 


0.000324 


-0.003488 


0.145197 


0.000535 



TABLE IV: Line strengths (a.u.) calculated in length L and 
velocity V forms for transitions between the selected odd- 
parity core-excited states with J = 1/2 and the ground 4s 
and singly excited 4ci3/2 states in W*^^ ion. 



Upper level 



First order 
L V 



RMBPT 
L V 



[4lj'll'j'{Ji)3l"j" (1/2)- 4si/2] transitions 

4si/24pi/2(l)3d3/2 0.002056 0.001818 0.002152 0.002151 

4p3/24d3/2(2)3d5/2 0.019622 0.017609 0.017775 0.017794 

4^1/24/5/2(2)3^3/2 0.004205 0.003780 0.003862 0.003865 

4p3/24d3/2(2)3d3/2 0.071111 0.063933 0.065253 0.065476 

4si/24d5/2(2)3p3/2 0.009202 0.008321 0.009124 0.009165 

4p3/24p3/2(2)3p3/2 0-001041 0.000942 0.001092 0.001093 

[40'4r/(Ji)3/"j" (1/2)- 4d3/2] transitions ~ 

4si/24d3/2(2)3p3/2 0.003858 0.003437 0.004065 0.004078 

4d3/24/7/2(3)3d5/2 0.006774 0.006079 0.006224 0.006230 

4d5/24/5/2(3)3d5/2 0.014631 0.013129 0.013505 0.013423 

4pi/24pi/2(0)3pi/2 0.003732 0.003355 0.003401 0.003407 

4d5/24/5/2(l)3d3/2 0.057817 0.051983 0.052786 0.052932 

4d5/24/7/2(2)3d3/2 0.006813 0.006127 0.006215 0.006210 



10" 
10" 
10'^ 
10" 

'cn 

10'° 

m 

2 10' -J 
.0 10° i 
10 
10= 
10= 

10' 



CO 
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6=4p4d('P)3d =P„^ - 4s 's' 
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Nuclear charge Z 



FIG. 1: Weighted transition rates between core-excited odd- 
parity states with J = 1/2 and 4s ^51/2 states as function of 
Z in Cu-like ions. 



follows the pattern of the corresponding calculation for 
Na-like ions in Ref. 2]. These matrix elements are cal- 
culated in both length and velocity gauges. 

Table IIIII lists values of uncoupled first-and second- 
order dipole matrix elements Z'^^\ Z^^\ B^^\ together 
with derivative terms p^^^"") for Cu-like tungsten, Z = 
74. For simplicity, only the values for the selected dipole 
transitions between odd-parity states with J — 1/2 and 
the ground 4s and excited ^3/2 states are presented. 
The more comprehensive set of data is given in Table IV 
of ^3 ) . The derivative terms shown in Table IIIII arise 
because transition amplitudes depend on energy, and the 
transition energy changes order-by-order in RMBPT cal- 
culations. Both length (L) and velocity {V) forms are 
given for the matrix elements. We find that the first- 
order matrix elements Z^'^ and Zy"^ differ by 10%; the 



L - V differences between second-order matrix elements 
are much larger for some transitions. The term p(dci'v) 
in length form almost equals Z^-^^^ in length form but in 
velocity form is smaller than Z*^^^ in length form by three 
to four orders of magnitude. 

Values of line strengths calculated as a square of cou- 
pled reduced matrix element 2] in length and velocity 
forms are given in Table HVl for the selected dipole tran- 
sitions between odd-parity states with J — 1/2 and the 
ground 4s and excited 4^3/2 states. A more complete 
set of data is given Table V of 0. Although we use 
an intermediate-coupling scheme, it is nevertheless con- 
venient to label the physical states using the jj scheme. 
We see that L and V forms of the coupled matrix ele- 
ments in Table HVl differ only in the third or fourth dig- 
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TABLE V: Wavelengths (A in A) and transition rates (Ar in s~^) for transitions from core-excited states QJ {Q = 
4:lj4:l' j' {J\)3l" j") to the the ground 4s and singly excited 4pi/2 states in W*^^ ion. Comparison with theoretical data obtained 
by using COWAN code from Ref. |20| . Numbers in brackets represent powers of 10. 



jj couphng 


Amcdf 


Armbpt 


^RMBPT 


^cowan 


j^cowan 


LS coupling 






transitions 








4p3/24d3/2(l)3d5/2 


5.919 


5.9248 


2.763[14] 


5.9110 


3.020[14] 


4s4/3d(3D) '*D3/2 
4s4/3d(373) 4i3i/2 
4p4d3d(3p) "^3/2 
4p4d3d(3p) *P3/2 
AsAfM('D) 2Pi/2 
4s4/3d(3D) 2p3/2 
4p4d3d(3p) 2p3/2 
4s4d3p(ip) 2p3/2 
4s4d3p(^P) 2Pi/2 


4j53/24d3/2(2)3d5/2 


5.912 


5.9200 


1.735[14] 


5.9069 


1.786[14] 


4Sl/24/7/2(3)3d3/2 


5.752 


5.7589 


1.863[14] 


5.7578 


1.832[14] 


4p3/24d3/2 (0)3^3/2 


5.744 


5.7512 


2.824[14] 


5.7501 


3.295[14] 


4P3/24d3/2 (2)3^3/2 


5.718 


5.7291 


7.Q26[14] 


5.7274 


7.954[14] 


4P3/24d3/2(2)3d3/2 


5.715 


5.7239 


8.195[14] 


5.7243 


9.738[14] 


4Pl/24pi/2(0)3p3/2 


5.615 


5.6347 


9.915[13] 


5.6419 


l.Q03[14] 


4Si/24d5/2(2)3p3/2 


5.230 


5.2409 


3.212[14] 


5.2322 


3.023[14] 


4Si/24d5/2(2)3p3/2 


5.222 


5.2328 


1.289[14] 


5.2252 


1.288[14] 




[4/i4;'j'(Ji)3/"j"-4pi/2] 


transitions 






4p4/3d(»P) *P3/2 
4p4/3d(3F) *Di/2 
4d4d3d(iG) ^D3/2 
4d4d3d(^P) ^Dg/a 
4p4/3d(ip») 2^3/2 
4p4/3d(3p) 2pi/2 
4p4/3d(iZ?) 2p,/2 
4p4d3p(3S) "1)3/2 


4Si/24p3/2(l)3p3/2 


5.912 


5.9175 


4.245[14] 


5.9093 


5.406[14] 


4Sl/24p3/2(l)3p3/2 


5.910 


5.9179 


2.314[14] 


5.9074 


2.855[14] 


4pi/24/7/2 (3)3^3/2 


5.752 


5.7544 


5.714[13] 


5.7407 


1.497[14] 


4^3/24^5/2 (l)3d5/2 


5.732 


5.7349 


1.680[14] 


5.7206 


6.059[14] 


4d3/24d5/2 (2)34/2 


5.723 


5.7325 


1.235[15] 


5.7256 


8.290[14] 


4d3/24d5/2(2)3d5/2 


5.722 


5.7267 


4.256[14] 


5.7210 


7.570[14] 


4d3/24d5/2(3)3d5/2 


5.719 


5.7268 


2.074[14] 


5.7130 


4.057[13] 


4pi/24d5/2(2)3p3/2 


5.233 


5.2431 


1.086[14] 


5.2354 


7.306[13] 


4pi/24d5/2(3)3p3/2 


5.230 


5.2400 


2.323[14] 


5.2548 


2.845[14] 


4p4d3p(3D) 
4p4d3p(='D) "Pi/2 


4pi/24d5/2(2)3p3/2 


5.229 


5.2395 


1.825[14] 


5.2546 


1.68Q[14] 
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FIG. 2: Weighted transition rates between core-excited even- 
parity states with J — 1/2 and 4p "P3/2 states as function of 
Z in Cu-like ions. 



its. These L-V differences arise because our RMBPT 
calculations start with a non-local Dirac-Fock (DF) po- 
tential. If we were to replace the DF potential by a lo- 
cal potential, the differences would disappear completely. 
The first two columns in Table HVl show L and V values 



of line strengths calculated in the first-order approxima- 
tion which differ by factor of 10. Thr last two columns 
indicate that including of the second-order contribution 
almost removes the L — V differences. 



III. RESULTS AND COMPARISONS WITH 
OTHER THEORY AND EXPERIMENT 

We calculate energies of core-excited even-parity 
4lj4l'j'{Ji)3l"j" as well as the singly excited 4/j states 
in Cu-like ions with Z = 30 -100. Reduced matrix ele- 
ments, line strengths, oscillator strengths, and transition 
rates are also determined for electric dipole transitions 
between the above mentioned core-excited and singly ex- 
cited states in Cu-like ions for the same range of Z. 

Table IVI lists theoretical data for selected transitions 
from core-excited levels to the ground 4s and a singly ex- 
cited 4pi/2 levels with largest values of transition rates in 
W''^+. Also, the comparison of theoretical data produced 
by the different methods and codes is inlcuded. Specif- 
ically, three values for wavelengths: Amcdf, Armbpt, 
and Xcowan are compared. The first values, Amcdf, were 
obtained as the first-order RMBPT values. It can be 
seen from Table IVl that the results obtained by COWAN 
code [i^l, Xcowam better agree with the MCDF results, 
Amcdf than with the RMBPT results, Armbpt- 

COWAN code j23| gives results which are generally in 
good agreement with experimental energies by scaling 
the electrostatic Slater parameters to include the corre- 
lation effects (Refs. iJl;]) and here the scaling factor of 
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TABLE VI: Wavelengths (A in A) and transition rates {Ar 
in s~^) for transitions from core-excited states QJ [Q — 
4lj4l' j' {Ji)3l" j") to the the ground 4s and singly excited 
Alj states in Cu-like Ce, W, and An ions. Comparison of 
theoretical data obtained by using RMBPT code with exper- 
imental measurements (Acxpt) from Ref. (a), Ref. Q (b), 
and Ref. ^ (c)). Numbers in brackets represent powers of 10. 



Transition Armbpt Af^^^'^ Acxpt 





Cu-Uke 


Ce, 


z58 












4si /t4 /T 


f3)3pi /o(5/2) 


4nQ /o 


9 


4558 


4 


16 


131 


9.452" 


4si /o4p-i/2 


(2)3si/2(5/2) 


4d3/2 


9 


9569 


8 


32 


12] 


9.940" 




(2)3p3/2(5/2) 


4d5/2 


10 


0518 


1 


58 


13] 


10.051° 


4j33/24d3/2 


(2)3p3/2(l/2) 


4p3/2 


10 


1114 


1 


96 


13] 


10.111" 


4d3/24d5/2 


(4)3p3/2(5/2) 


4d3/2 


10 


1435 


3 


52 


13] 


10.141" 


4p3/24/7/2 


(2)3d3/2(3/2) 


4pi/2 


11 


3223 


8 


62 


11] 


11.321" 


4Pl/24p3/2 


(2)3pi/2(5/2) 


4d3/2 


11 


4202 


7 


77 


12] 


11.412" 


4d3/24d5/2 


(l)3d3/2(l/2) 


4pi/2 


11 


4912 


1 


24 


14] 


11.491" 


4d5/24/5/2 


(3)3d5/2(3/2) 


4d3/2 


11 


8120 


1 


64 


13] 


11.816" 




Cu-like 


w, z= 


74 












4P3/24d5/2 


(3)3d3/2(5/2) 




6 


8197 


2 


83 


12] 


6.816'' 


4p3/24/7/2 


(3)3d3/2(5/2) 


4/7/2 


6 


8269 


4 


26 


12] 


6.827*' 


4Pl/24p3/2 


(2)3d3/2(l/2) 


4pi/2 


6 


8577 


2 


32 


12] 


6.858*' 


4Si/24p3/2 


(l)3d3/2(3/2) 


4Sl/2 


6 


8885 


3 


93 


12] 


6.884*' 


4d3/24d5/2 


(4)3d3/2(5/2) 


4/7/2 


6 


8947 


1 


89 


12] 


6.896'' 


4p3/24d5/2 


[3)3d5/2(7/2) 


4d5/2 


7 


0724 


3 


41 


13] 


7.075'' 


4P3/24/7/2 


(3)3d5/2(5/2) 


4/7/2 


7 


1095 


1 


31 


12] 


7.108'' 


4»3/o4f7/o 


3)3d5/2(9/2) 


4/7/2 


7 


1117 


6 


78 


12] 

J 


7.113'' 


4r>Q /o4 fr. /o 
1 2.^J of Z 


f3~)3ds/of5/2") 


7/2 


7 


1302 


1 


96 


13] 


7.131'' 




^-LyoLt5/2V ' / ^1 


4dr /r, 

^1^5/2 


7 


1354 


1 


59 


131 


7.137'' 


^^1/2^/7/2 


^'±jO"5/2lO/ 


™5/2 


7 




1 

L 


fii 

Ul 


191 


7 OztO*" 




l^joa5/2(,0/ A) 


4P3/2 


/ 


0/1 70 



Z 


UD 




7 9/1 8^ 


/I ^ /If 
4Sl/24/7/2 


3j3d5/2(3/2j 


/I ^ 

4d5/2 


I 


zOUo 


1 
1 


00 


13] 


1 .ZoZ 


4pi/24d3/2 


(l)3d3/2(5/2) 


4d3/2 


7 


2916 


1 


64 


13] 


7.293'' 




Cu-like 


Au, 2'= 


-79 












4d5/24d5/2 


(4)3p3/2(7/2) 


4d5/2 


4 


3974 


1 


65 


14] 


4.39^= 


4d3/24/7/2 


(4)3p3/2(9/2) 


4/7/2 


4 


4690 


1 


79 


14] 


4.69" 


4d5/24/7/2 


(2)3d3/2(5/2) 


4d3/2 


4 


7690 


1 


66 


14] 


4.76" 


4/7/24/7/2 


(4)3d3/2(7/2) 


4/7/2 


4 


7997 


4 


53 


15] 


4.80" 


4/5/24/7/2 


(4)3d3/2(7/2) 


4/5/2 


4 


8214 


1 


98 


15] 


4.82" 


4d5/24d5/2 


(4)3d5/2(5/2) 


4p3/2 


4 


9373 


2 


51 


14] 


4.93" 


4d5/24/5/2 


(3)3d5/2(3/2) 


4d3/2 


4 


9610 


2 


58 


14] 


4.96" 


4Sl/24p3/2 


(2)3p3/2(l/2) 


4P3/2 


5 


2208 


1 


15 


13] 


5.22" 


4d3/24d5/2 


(3)3d3/2(7/2) 


4/5/2 


5 


7207 


2 


08 


12] 


5.72" 


4Si/24p3/2 


(2)3d5/2(3/2) 


4Sl/2 


5 


9029 


3 


93 


13] 


5.90" 


4P3/24d5/2 


(l)3d5/2(7/2) 


4d5/2 


5 


9427 


2 


10 


13] 


5.94" 


4p3/24d3/2 


(2)3d5/2(7/2) 


4d5/2 


6 


0107 


1 


28 


12] 


6.01" 


4Sl/24d5/2 


(3)3d5/2(l/2) 


4p3/2 


6 


0801 


1 


61 


13] 


6.08" 



0.85 was used. Also, in COWAN code the LS coupling 
scheme is implemented and corresponding LS designa- 
tions were added to jj designations in all tables which 
include comparisons with theoretical data produced by 
COWAN code. 

A comprehensive set of theoretical data and compar- 
isons for transitions with largest values of transition rates 
in W*^"*" similar to presented in Table is given in Ta- 
ble VI of the accompanying EPAPS document [H. It 
includes transitions from core-excited levels not only to 
the ground 4s and to one singly excited level 4pi/2 as in 



Table Ivl but to all singly excited Alj levels. 

The similar comprehensive set of theoretical data and 
comparisons for transitions between core-excited states 
4lj4:l'j'{Ji)3l" J to the ground 4s and singly excited 4L 
states is presented for Ce^^+ ion in Table VII of Ref. [Ig 
and for Au^"+ ion in Table VIII of Ref. ^3 ■ Comparisons 
with theoretical data obtained by using RMBPT, MCDF, 
and COWAN codes show that the difference in results can 
be explained by the second-order corrections to energy 
and dipole matrix elements included in RMBPT. 

The trends of the Z-dependence of transition rates for 
the transitions from core-excited odd-parity states with J 
= 1/2 to the ground state 4s states are presented in 

Fig. n . The trends for the transitions from core-excited 
even-parity states with J — 1/2 to the singly excited state 
Ap ^^3/2 are shown in Fig.|2| More transition rate figures 
are given in [I8] . 

We see from the graphs that transitions with smooth 
Z-dependence are rarer than transitions with sharp fea- 
tures. Smooth Z-dependence occurs for transitions from 
doublet and quartet core-excited states. Usually, singu- 
larities happen in the intermediate interval of Z = 40 - 
60 when neither LS nor jj coupling schemes describe the 
states of these ions properly. One general conclusion that 
can be derived from these graphs is that the smooth Z- 
dependence occurs more frequently for transitions from 
low-lying core-excited states. 

Singularities in the transition-rate curves have two 
distinct origins: avoided level crossings and zeros in 
dipole matrix elements. Avoided level crossings result 
in changes of the dominant configuration of a state at a 
particular value of Z and lead to abrupt changes in the 
transition rate curves when the partial rates associated 
with the dominant configurations below and above the 
crossing point are significantly different. Zeros in tran- 
sition matrix elements lead to cusp-like minima in the 
transition rate curves. Examples of each of these two 
singularity types can be seen in Figs. ^ and [3 as well as 
in Figs. 1-10 presented in EPAPS document 18]. 

In Table wavelengths and electric-dipole transition 
rates along with comparison with experimental data are 
presented for transitions in Cu-like Ce, W, and Au. The 
table is limited to identification of experimentally mea- 
sured transitions given in Refs. Sid. A more com- 
prehensive set of theoretical data assigned to each ex- 
perimental transition is presented in Table IX of EPAPS 
document 0. We mentioned previously that all possi- 
ble 4/iji4?2j2(>/i)3/3j3(J) - 4:lj transitions produce 2294 
spectrum lines. These lines in Ce^^+, W'^^, and Au^""*" 
cover the spectral regions from 6.8 - 21.9 A , 3.5 - 9.8 
A , and 2.9 - 8.7 A , respectively. A number of spectral 
lines becomes smaller by a factor of ten when we con- 
sider transitions with largest values of transition rates, 
gAr- The number of transitions with gAr > 10^^ s^^ is 
about 200 for Ce^^^, and a number of transitions with 
gAr > lO^-* s-i is about 130 for W^^+ and about 190 for 
Au^'^+. Those transitions with the largest values of gAr 
cover the spectral regions smaller by factor of 3-4 than 
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FIG. 3: Synthetic spectra of Cu-like W calculated using 
the RMBPT atomic data and Gaussian profiles with (AA 
= 0.02 Angstrom). The scale in the ordinate is in units of 



These atomic data are important in modeling of M- 
shell radiation spectra of heavy ions generated in elec- 
tron beam ion trap experiments and in M-shell diagnos- 
tics of plasmas. For example, x-ray M-shell spectra of 
W ions collected at the different energies of the electron 
beam at LLNL electron beam ion trap indicate the ex- 
istence of strong Cu-like spectral features 'Q^. Also, the 
distinct features of x-ray M-shell spectra of W ions in 
a spectral region from 5 to 6 A produced by labora- 
tory plasmas are the brightest Ni-like lines and Co- and 
Cu-like spectral features. In particular, Cu-like satellite 
lines constitute most of the Cu-like spectral features, and 
Cu-like autoionization levels make an im por tant contri- 
bution in ionization balance calculations |23 |. Synthetic 
x-ray spectra of Cu-like W computed with different res- 
olution of 0.02 A and 0.002 A are presented in Fig. |31 
and 0] respectively. It was assumed that spectral lines 
have the intensities proportional to the weighted transi- 
tion probabilities and are fitted with the Gaussian profile. 
Specifically, the spectrum in Fig. |21 includes 3^ — il' tran- 
sitions and covers the spectral region from 4 A to 6.4 A. 
The most intense peaks at about 5.72 A and 5.9 A are 
formed by 3d — 4/ transitions and are shown with better 
resolution in Fig. 0] 



300 
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FIG. 4: Synthetic spectra of Cu-like W calculated using 
the RMBPT atomic data and Gaussian profiles with (AA — 
0.002 Angstrom). The scale in the ordinate is in units of 



that mentioned for all lines. We can see that the number 
of predicted spectral lines even with the largest values 
of gAr still much larger than a number of experimental 
lines shown in Table [711 Also, the interval between these 
lines is about 0.01 A or less which is comparable with the 
accuracy of experimental measurements. In this case, it 
could be reasonable to assign not a single transition as in 
Table IVII but 3-5 transitions to identify the experimen- 
tal peak, as it is demonstrated in Table IX of EPAPS 
document M. 



IV. CONCLUSION 

In summary, a systematic second-order RMBPT 
study of the energies and transition rates for 
[4:liji4:l2j2{Ji)3l3j3 (J) - 4:1 j] electric-dlpolc transitions 
in Cu-like ions with the nuclear charges Z ranging from 
30 to 100 has been presented. The retarded El matrix 
elements included correlation corrections from Coulomb 
and Breit interactions. Both length and velocity forms 
of the matrix elements were evaluated and small differ- 
ences (0.4 - 1 %), caused by the non locality of the start- 
ing DF potential, were found between the two forms. 
Second-order RMBPT transition energies were used in 
the evaluation of transition rates. These calculations 
were compared with other calculations and with avail- 
able experimental data. For Z > 40, we believe that the 
present theoretical data are more accurate than other 
theoretical or experimental data for transitions between 
the 4^iji4Z2j2(^i)3^3j3 (J) core-excited states and the 4lj 
singly excited states in Cu-like ions. The results could 
be further improved by including third-order correlation 
corrections. These calculations are presented as a theo- 
retical benchmark for comparison with experiment and 
theory. In addition, the application of generated atomic 
data to modeling and interpreting of x-ray M-shell spec- 
tra of heavy ions is discussed. 
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